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Abstract 

A hybrid model for baryons based on a dynamical interplay between rel- 

ativistic three-quark bound states and soliton configurations of mesons is 

constructed. The Bethe-Salpeter equation for diquarks and the Faddeev 

equation for diquark-quark bound states in the background of a soliton is 

solved. The results show that baryons are very much like hybrids containing 

both, solitonic meson clouds and three-quark correlations. 
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I. INTRODUCTION 



At present there exist two generic classes of models to describe baryons. On the one 
hand, there is the picture of baryons as chiral solitons. The soliton picture is based on 
considering QCD for an arbitrary number of color degrees of freedom, Nc- In the infinite 
color limit, Nc oo, QCD reduces to an effective theory of infinitely many weakly inter- 
acting mesons 0. Although this effect ive meson theory cannot be constructed explicitly, 
Witten conjectured that within this theory baryons emerge as solitons 0. Based on this 
conjecture phenomenological effective meson theories have been developed which possess 
soliton solutions. The most prominent is perhaps the Skyrme model In the limit of 

an infinite number of colors the soliton description is the only model for baryons. On the 
other hand, for a finite number of colors a baryon is customarily considered as a bound 
state of three valence quarks. Such valence quark models are motivated by high energy 
scattering experiments which have revealed a partonic substructure. Starting with these 
experimental facts many models, which are based on the valence quark picture, have been 
developed. These models include the non-relativistic quark models and its relativis- 
tic extensions 0, parton models which are directly based on the scale invariance [|r^, bag 
models and diquark-quark models []T2|,|TB[. 

The valence quark picture directly leads to the quantum numbers of a physical baryon 
whereas the soliton can only be interpreted as a baryon with good spin and flavor quantum 
numbers after collective quantization. Yet, baryons as solitons are conceptually better 
suited for the description of low energy properties because they straightforwardly embed 
the useful feature of chiral symmetry and its spontaneous breaking. In any event, despite 
their successes both pictures possess only limited ranges of applicability. 

Since the advantages of both pictures are in some sense complementary a unification of 
the two approaches seems desirable. In principle, the chiral bag model [0,0 represents 
such a combination since inside the bag it contains explicit quark degrees of freedom 
whereas a chiral soliton field surrounds the bag. As a consequence of the Chesire cat 
principle [^| experimental measurable quantities should not depend on the radius chosen 
for the bag. Explicit calculations show that the Chesire cat principle is not universally 
valid. In particular, the singlet axial matrix element depends strongly on the bag radius 
|l^,18|. Therefore, a hybrid model is desirable which connects both classes of models in a 
dynamical fashion. In this context the Nambu-Jona-Lasinio (NJL) model [19| is for the 
moment unique. On the one hand, it possesses soliton solutions pO| of meson fields which 



themselves are obtained as bound antiquark-quark states |21|. On the other hand, within 
this model baryons may be described as three-quark bound states via the use of diquarks 
22-p8[, while the meson fields are fixed to their vacuum expectation values. In particular 



this model is unique because with the help of path integral hadronization techniques 
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a consistent unification of botli approaclies is possible witliout any double counting of 
correlations. 

Because of the enormous computational effort needed, the realization of a rigorous 
self-consistent solution is for the moment not feasible. Nevertheless, an approximate 
evaluation of this hybrid baryon can be accomplished within a four step procedure to be 
carried out in this paper: In the proceeding section the transformation of the NJL model 
with a pointlike interaction of color octet current into an effective theory of mesons, 
diquarks and baryons is described. For completeness we briefly repeat the first two steps 
in the beginning of the third section: First, we construct a static ground state solution in 
the absence of diquark and elementary baryon fields. Second, we solve the Bethe-Salpeter 



equation for a scalar diquark in the solitonic background field ||3y]. In addition, we derive 



the Faddeev equation for arbitrary quark propagators. In section |I^ we discuss numerical 
results for the solutions of the Bethe-Salpeter equation as well as of the Faddeev equation. 
In section we employ these results to define the hybrid model and discuss some static 
properties. We close with conclusions and a outlook in section |V|. Some details of the 
calculation and a few lengthy formulas are left to appendices. 



II. HADRONIZATION OF THE NJL MODEL 

As stated in the introduction we consider a NJL model for two flavors 

CNJL = q(i^-m')q-\gf,fa (2.1) 

with a pointlike interaction of color octet flavor singlet currents = gA^7^g/2. Here 
q denotes the quark spinors, rhP = diag(mQ, mg) the current quark mass matrix for two 
flavors and Xci^ ~ ••' ~ 1) generators of color SU{Nc)- The interaction part 

can be Fierz-rearranged to solely attractive channels 

Ant = ^1 (l - ^) (qAaq) (gA"g) + ^ (gF^g^) (gT"g) , (2.2) 

where q"^ = Cq^ denotes the charge conjugated Dirac spinor. Furthermore, we have 
defined the vertices 

2 ^ ' ^ \V2j^ 2 



A° = l,— O^, r"=^ — with a = {A,a,a) (2.3) 



for the quark- ant iquark and the quark-quark interaction, respectively. corresponds to 
the set of Dirac matrices 

c III jYl5i 
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and t'^ are isospin matrices. Because of the Fierz-transformation the couphng constants 
are restricted to gi = g2 = g, which is necessary to preserve consistent Nq counting. 
Note, that the diquark channel is suppressed by a factor I/Nq in comparison to the 
meson channel. Henceforth we will confine the discussion to the physical case Nc = 3, 
unless explicit noted. 

To convert the pure quark NJL model ( |2.1D into an effective hadron theory^ we 
introduce collective meson and baryon fields and into the generating functional 



Z = J VqVqexp (^j Cnjl 



via the identities 



j l^xJixa - qKq) = J Pxa'E'v^^exp (^i J ip'^ixa - qKq) 

I VBiVBi5 - ^(gT„g)P^) 5 (pf - ^-fqP^qV^q 



VBEVB^V'^'^V'^'^ 



a a 



(2.5) 



(2.6) 



(2.7) 



exp t 



B^-'^-f (gT„g) P^q] + [b^ - ^qP^ (qV^q^) 



Note, that we are working in Euclidean space where we have used the Wick rotation 
t — > — ir. For notational simplification we use the abbreviation J = J di^x = J d^rdr in 
the exponents. The operator P^ projects the three quarks onto the quantum numbers of 
the considered physical baryon states. In particular, the color of the third quark has to 
be chosen to build a colorless baryon wavefunction. 

As intermediate building blocks for the baryon field we introduce also diquark fields 

1 = j Vk^Vk*J{k*^ - gr„g^)5(«:„ - gT.g) (2.8) 

= J I?/€„P<PA°PA"*exp Q J [« - gr„g^)A" + A"*(«:„ - gT.g)]) . 



Replacing all terms of fourth order in the quark fields with help of the constraints (2.6) 
and ( [^.8| ), and integrating out the auxiliary fields Xa, f^a and , one obtains the following 
form for the generating functional 



VqVqD'^V^VM)BVBVm''Vm' 



(2.9) 



exp 



q { i-i^d^ - m° - + \ q + gA^^" + ^°A„g 



exp j ^gAg^ + Ig^Ag + ^tr^^ ^ AtrAA - i + B^^mf^ j 



^Similar hadronization approaches have been considered in pl|j3 
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where we have introduced the compact matrix notation for the meson and diquark fields 
V9 = v^aA", A = r„A", A = r„A"* as well as for the baryon sources = ^^P^ and 
_ p/3i|ra_ xhe symbol tr corresponds to the trace over color, fiavor and Dirac spinor 
degrees of freedom. 

To eliminate the quark degrees of freedom we are working within the Nambu-Gorkov 
formalism developed originally in the theory of superconductivity P5|-p5[|. For that pur- 
pose we introduce combined Grassmann fields^ 



q 



for the quark fields and 



[2.10) 



J 



[2.11) 



for the baryon sources. Now, we integrate out the quark fields q and q with help of the 
Nambu-Gorkov formula 

VqVqexp {qg-\ - ^qAq' - ^(fAq + [gA*°^„ + ^aA"q 



PqPqexp ( - qG ^q + q^ + ^q 



DetG-M'exp ( |^G^ 



[2.12) 



where we have defined 



/ 



v 



and the inverse quark Greens function 



[2.13) 



G 



-1 



V 



-CA* 



-AC 



g = v^g^v. 



[2.U) 



The off-diagonal elements of G are the so-called anomalous Greens functions which are 
related to the amplitude for adding or subtracting a pair of quarks to the system. The 
transformation operator V = JG, which we have introduced for technical reasons, is a 
combination of the self-adjoint unitary transformation J = 1/3'-/^ and the G-parity operator 
G = e'^'^^/'^G. 



^Quantities in the Nambu-Gorkov formalism are denoted by boldface letters. 
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The normal quark Greens function is represented by 

g-^ = i^d^ + = ^Q-i + l^i|r"^". (2.15) 

3 3 



From eq. (|2.13| ) we observe that the baryon sources are contracted with the diquark 



fields Aa- Hence these sources only couple to a single quark. Since in the ladder approxi- 
mation (described below) a three quark bound state cannot be aff'ected by such a couphng 
it is sufficient to restrict the quark Greens function to Qq^. 
Finally, the generating functional is given by 

exp {A[^, a, a, VI/"] + , (fr^E^ + 5f vl/^)) (2.16) 

with the effective action 

^[V?, A, A, = (2.17) 

A,[v, A, A, + Aa/b, A, A, + An.lv] + A4A, A]. 

The first part, the so-called quark determinant, 

= ^TrlogG-^ (2.18) 

carries the full information of the underlying quark spectrum via the functional trace Tr. 
The second part 

A^al = i*"(G^'^%;3*^ (2.19) 

contains the valence quark part of the baryon propagator (G™^)q,/3 originating from the 
inversion of the quark propagator Go (c/. appendix B for the derivation of (G^"')^/?.). 
Both residual terms 



tr(^2 



8gi J 

Ad = f d^x [trAAl (2.20) 



8^72 

are pure mass terms for the meson and diquark fields, respectively^. 



^Note the additional factor Nq = 3 in the diquark mass term Ad as compared to the mass term 
Am of the mesons. 
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III. TWO AND THREE QUARK CORRELATIONS IN A SOLITONIC 

BACKGROUND 



A. The static soliton background 



In this section we are interested in the behavior of quark-quark and diquark-quark 
bound states in the background field of a sohton configuration. In the first step we 
neglect the explicit baryon sources. With this restriction it is possible to expand the 
effective action ( |2.17|) 



^[V?, A, A, = 0, = 0] = [i^] + A^^^ [<^, A, A] + ^(2) A, A] + . . . 
up to second order in the diquark fields where (see also eqs. ( p.l9| ) and ( |2.2q )) 

^«[^,A,A]=4i), 
^(2)[y.,A,A]=4^)+^,. 



(3.1) 

(3.2) 
(3.3) 
(3.4) 



Within Schwinger's proper time regularization description ||36| the parts stemming from 
the quark determinant ( p.l8|) can be written as 



4 yi/A2 s 



A^P = -Tr / dsA^e-'^\ (3.5) 
^ 4 7i/A2 

Af^ = -Tr / ds A2e-'^° - -Tr / dss / rfaaAie"^"^" Aie^^^^-")^". 
^ 4 Ji/A2 4 7i/A2 Jo 



1/A2 4 J1/A2 

The Nambu-Gorkov matrices Aj originating from the expansion of the operator 
(G"^)^G"^ = Ao + Ai + A2 have the explicit form 



An 



Ai 



A, 



{g-'Vg-' 
v-^{g-^yg-^v 



cA^g-^ + v-\g~^yvcA 

( AtA ^ 
-CA^AC 



V 

/ 

V 



-{g-yAc - K^cv'^g-^v ^ 







(3.6) 



(3.7) 



(3.8) 



By expanding ^^^^ to quadratic order in we can relate ||2T| , |37|1 the parameters (5^1, 
mo = ttlq = mQ and A) of the model to the pion mass m^r = 135MeV and the pion decay 
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constant f-j^ = 93MeV, leaving undetermined just one parameter which we choose to be 
the constituent quark mass m. This mass is the vacuum expectation value of the scalar 
field and refiects the spontaneous breaking of chiral symmetry. In principle the diquark 
coupling constant g2 is fixed by the Fierz transformation, however, we will relax the 
restriction g2/gi = 1 thus effectively treating g2/gi as a parameter to study the influence 
of the diquarks. 

Note, that for each term in eq. ( |3.5| ) the Nambu-Gorkov trace can be worked out 
analytically. In addition, due to the diagonality of Aq and the vanishing trace of Ai it 
is obvious that A = A = is always a solution of the equation of motion though it is 
not necessarily the solution of least action. For the remaining functional trace we use the 
color degenerated eigenstates of the inverse propagator for quarks in the background of a 
static meson fleld (p 

g~'^\n, v) = {i<^ - !f) \n, v) = P {-dr - h@) \n, u) = (3 {iun - e^) |n, v) (3.9) 

and perform an additional color trace tr^. In analogy, the charge conjugated quark 
propagator Qq^ can be worked out similarly: 

v) = - ifY V\n, v) = (3 {dr - he) \n, v) = (3 {-iun - e^) |n, v). (3.10) 

For the meson flelds we make use of the polar decomposition 

^ = ^1]""' =7X1^, (3.11) 

where we have flxed the chiral radius $ to its vacuum expectation value m. Additionally, 
we adopt the well-known hedgehog ansatz for the chiral fleld 

f/ = exp(iT-fe(r)). (3.12) 

Due to the static fleld conflguration it is possible to extract from the zero order term 
A^^^ = —TEsoi[Q] (|3.2| ) an energy functional 

Esol[Q] = Esea + E^al (3.13) 
3 foo 1 „ r 

Esea = - 1^^ ds^== E e-''^ + <fl / - ^ose{r)) 



3 

E^al = -^^yaiil + sgn(e„a/)). 

Neglecting quark-quark correlations, the valence quark contribution E^ai is composed of 
the color degenerated single particle energies of the valence quark level. This state is 
the only one which is bound by the soliton, i.e. —m < e^ai < itl. The soliton conflgura- 
tion is characterized by the chiral angle Qsci^) which is self-consistently determined by 
extremizing the energy functional -^[O] []38| - 
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B. Bethe-Salpeter equation for diquarks 



In the following we derive the equation of motion for the diquark fields in the solitonic 
meson background field. The solutions of this equation define the eigenmodes 6^ and 6* 
of the diquark fields 



A = b,,{A{r,t)TCV-'r = b^r%^{r,t), 
A = {VCTA{r,t)rb: = b*J",^^{r,t) 



(3.14) 



which we have redefined including the matrices V and C into the diquark vertices. For 
the present qualitative discussion it is sufficient to consider only an S-wave scalar diquark 
field. In respect of the Pauli principle the only possible ansatz for the vertex functions 
^Ur,t) is 



r--|f.7a (3.15) 



with Xc{a = 2, 5, 7) being the antisymmetric Gell-Mann matrices of the color group. 
Using the Fourier-transform in the time coordinate 

duj 



A(r,T) 



27r 



A(r, iuj)e 



(3.16) 



and the spatial representation \E'jy(r) = (r|z/) for the quark eigcnstatcs the residual trace 
can be worked out. The second order term of the effective action, A^g'\ together with Ad 
determines the Bcthc-Salpeter equation as the equation of motion for the diquark fields 
SA^^ySA*{r, -icu) = 0: 



0. 



[y dr'r'^D^^^{r,r';iu)A^{r',tu) 
The inverse diquark propagator 

Dliir, r'; iu;) = K{r, r'; icv) - i^Kr - r') 
is expressed in terms of a local mass term and a bilocal kernel 

K{r',r-iu) = Y.j dQdn'Riiu; 6^,6^) [*t (r.')*^(r')*i(r)^r.(r)' 

where 

R{iuj; e^, 



(3.17) 



1 e 



2e. 



16^/^ 

(up' + (e^ - e^) 



2' A2 



2 e„ 



2' A2 



(3.18) 



(3.19) 



(3.20) 



1 F (|, [ae^ + (1 - a)e^ + a(l - a)cjVA2 
daoL — — 



ae^ + (1 — Q;)e2 + a{\ — a)u!'' 



3/2 
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in the proper-time regularized version of the quark loop. 

For the ongoing discussion a normalization of the diquark field Aa{r, iu) is necessary. 
In principle such a normalization would be obtained in the framework of second quan- 
tization. Equivalently, we demand that the total baryon charge, to which the diquarks 
contribute, equals unity. For that purpose we first calculate the expectation value of the 
baryon number operator B = Ic/Nc = lc/3 (see appendix C for details.). In analogy to 
the effective action ( p. 171 ) the expression ( |C5| ) is expanded in terms of diquark fields. The 
leading order yields the value 1/3, which is nothing but the contribution of a single quark. 
Therefore, we require for the second order term the normalization condition 

Bdig = I (3.21) 

with Bdig defined in ( |C1CI| ). As the solutions of the Bethe-Salpeter equation appear al- 
ways in pairs ±iU(iiq, the normalization condition ( p.21| ) allows us to distinguish between 
diquark and antidiquark solutions. 



C. Faddeev equation for a quark-diquark bound state 



Up to this point, we have dealt with the generating functional ( p. 16 ) for meson, diquark 



and baryon fields. Below we derive an effective meson-baryon action generalizing the 
approach of P^. For this purpose we perform the diquark integration. Since the diquark 



fields appear to all orders we introduce external diquark sources j,j and rewrite Z ~ 
JVy^V^'^V'^'^Zdigij,]] (^J^) in the form 

Zdic[j, j] ~ / VAVA exp (^(2) ^] _^ A^nt[A, A] + z [Aj + jA] ) 

~ exp ^An^Jj, J-]) / 1? ADA exp (AD.-JA + ^(Aj + jA)) . (3.22) 
The integration over the diquark fields results in 

Zdiq[j,j] ~ exp ^Aa/[i, T^rlj exp [jDdiqj) . (3.23) 

Note, that the interaction part Amt = -Aval is completely determined by the valence quark 
part of the baryon propagator. Using the geometric series 



(i + ^oA^~oA)"' = (i + Gob^ri^gA^bi] 

oo 



a 



a=0 

to reexpress Gn in eq. (P^ we obtain the baryon propagator 
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(a = 0) 



+ 



+ 



(a = l) 



+ 



+ 



(a = 2) 



+ 




+ 



+ 



+ 



+ 



(a > 3) 



FIG. 1. Feynman diagrams of the series (|3.25| ) for the baryon propagator G^'^ 



y-'B )al3 



GoDaS^digGo^diaD-ypSo + • • • • 



(3.25) 



as an infinite series in quark Qo and diquark propagators Dajs- This series is drawn as a 
sum of Feynman diagrams in Fig. ||. For the calculation of bound states it is only necessary 



to consider the quark exchange graphs in the series ( 3.25 ). Neglecting the self-interaction 
graphs is justified because we have an effective interaction and therefore all self-interaction 
graphs can be absorbed in the effective coupling constants. This simplification is nothing 
else than the well-known ladder approximation which has the advantage that the baryon 
propagator can summed up completely: 
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dig^OJ- diqD'l-n^diqGo^'diq^€0't 

l-{Gl)^sH^'y\Gl),p. (3.26) 



Here we have introduced the free quark-diquark propagator 

{G%)al3 = GoDal3 (3.27) 

and the quark exchange operator 

H^' = Tl^GoTi^. (3.28) 

We are now able to integrate over the baryon source fields \E'",'i'". This leads to an 
effective baryon action 

^?Jry = B-iGl^'rjBf" (3.29) 
for the baryon fields B°' = P'^B^ und -B" = B'^P'^ . The inversion of the baryon propagator 

(<^B^)q/3 = (<^ij)a/3 - Hal3 (3.30) 



follows directly from eq. ( |3.26D . Using 



- (3.31) 



55" 

one finally arrives at the Faddeev equation in operator form: 

{GtUB^ = 0. (3.32) 

In the next step we construct the kernel of the matrix {G'"^'')'^^. For this purpose we 
expand the baryon field 

B{x,y)= dujdEY^ ^ aa,,(E, u;)5°^(x, X4; 1/4) (3.33) 
in eigenfunctions 

Bl{x, X,; y, y,) = J—A^{x)e-'^^^^.{y)e-'^y- (3.34) 
(27r)" 

of the free quark-diquark propagator 

{GD-^Blix, X4; y, 2/4) = -/3 {^E - e,) [u' + ul) B^x, x,- y, y^). (3.35) 
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The projection operator ( p.7|) on the third quark projects out the positive energy 
eigenstates {e^, > 0). Transforming the energy variables E, uo to the total energy = E+u 
and the relative energy Cl = ^{u — E) we obtain the Faddeev equation in terms of the 
coefficients aauii^ji^)'- 

■1 



{iQ, iQ) 



(3.36) 



+ 



4 





p/3 

diq 











0. 



The summation runs over the one-particle quark eigenstates \fi) and \k) of Hq. To reduce 
the numerical effort we restrict the diquark states to the energetically lowest state T^ig = 
Tdi^- In fact, this is the only bound state; all other states lie above the two quark 
threshold and, if at all, influence the quark-diquark bound state only weakly. 

It is important to equip the absolute and relative energies, Q and Q, respectively with a 
physical meaning. These quantities do not yet have apparent interpretations since due to 
the presence of the static and localized soliton Lorentz covariance is lost. Hence we are still 
lacking an energy scale. In order to deflne a physically relevant energy scale we consider 
the limit of vanishing diquark-quark coupling (F^jg —^0). In this limit the different quark 
states decouple and the Faddeev equation reduces to a purely algebraic equation. For the 
valence quark part we have the three solutions iQ = 2{iCl + e^ai) and iQ = 4(— zfi ± Udiq). 
When the energy transfer vanishes, i.e. Cl = 0, we should obtain the physically reasonable 
solution iQ = ojdiq + Cmi for the energy of a free baryon. This result can only be realized 
if we enforce the variable substitution — > 17 + |(e„a/ — ^diq)- Thus, two of the three 
solutions join the physical solution. The third one yielding a negative solution is rejected 
because we have projected to positive energy states. Since the limit of vanishing diquark- 
quark coupling should continuously emerge from the complete solution of the Faddeev 
equation it is necessary to perform the variable substitution^ 17 ^ + |(e^ai — Udiq) also 
in the Faddeev equation (|3.36|) . 

For the numerical solution of this integral equation^ we discretize the integration 
variable VLj = jACl; j = 0, . . . ,N, yielding the kernel as an x matrix. We then 
generalize the matrix equation to the eigenvalue equation 



^In the following we suppress all diquark indices, e.g. aa^f_i{i^,iCl) = a^(iil,iS7). 
similar variable substitution was implemented by Ishii et. al 

similar approach was employed for the solution of the Bethe-Salpeter equation for the 
diquarks pO| . 
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FIG. 2. Dependence of the soliton energy from the valence quark energy for three different 
parameterizations of the soliton profile. Results are for m = 450MeV. 



iql 



2 

dig 







\k){k 






+ ^^diq + — €yal 



afj^{iQ,iflj) = Xj{iQ)a,j{iQ,iQj) . (3.37) 



A solution of the Faddeev equation is given for a vanishing eigenvalue Xj{iQB) = with 
the energy Qb- Within this formulation there exists only one positive energy solution 
with Qb < ^vai + ^^diq- Possible negative solutions will again be rejected. If we replace 
Tdiq in ( |3.37D with x^diq we can explore the dependence of the solution on the coupling 
between quark and diquark. In this manner we observe that Qb approaches the free 
solution eyai + ^diq for decreasing x- This justifies the above described procedure. 



IV. NUMERICAL RESULTS AND DISCUSSION 

In the first step we calculate^ the chiral angle Qscij") self-consistently in a box with 
radius D by extremizing the energy functional ( |3.13| ). In order to later estimate the 
influence of the quark-quark and the diquark-quark bound states on the soliton we in- 
troduce a dimensionless parameter a, which measures the extension of the soliton profile, 
via 0(r) = Qsci'f^/a). This choice is motivated by the softness of the scaling mode, i.e. a 



'''A complete description of the numerical treatment is given in 
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0.0 759 934 175 

0.2 761 934 173 

0.4 758 929 171 

0.6 672 808 136 

0.8 465 560 95 

1.0 295 360 65 

1.2 161 208 47 

TABLE I. The energy w^jg of the diquark in the soliton, the energy 2€^ai of two uncorrelated 
quarks and the binding energy of the diquark Bqq = 2eyai — ^diq in dependence on the scaling 
parameter a for constituent quark mass m = 450MeV and g2 = Qi- 

weak dependence of the total soliton energy on the parameter a in comparison to other 
variables. The energy of the valence quark, on the other hand, depends sensitively on 
a. As shown in Fig. |^ the soliton energy for the scaling mode is almost constant in the 
displayed range. For comparison the dependencies for two alternative parameterizations 
are presented. Both profiles are chosen such that they not only fulfill the usual boundary 
conditions but also the derivatives at the boundaries are unchanged. 

In the next step we solve the Bethe-Salpeter equation (|3.17|) for the diquark fields^ In 



Tab. I we display the solution of the Bethe-Salpeter equation for different soliton profiles. 
The important result is that the diquarks are bound {ujdiq < '^^vai) for all considered 
profiles. Furthermore, the binding energy Bqq = 2evai—^diq decreases with growing soliton 
size. For a small soliton extension (a < 0.4) the diquark energy is almost independent 
of the soliton. Note, that the limit a — > corresponds to a diquark in an empty box. 
On the contrary, for larger soliton extension (a > 0.5) we observe a strong dependence 
on the scaling factor a. This behavior of the diquark energy is inherited from the one of 
the valence quark energy. This indicates that the properties of the diquark are essentially 
determined by the valence quarks. 

A further interesting point is the increasing diquark binding energy with growing con- 
stituent quark mass when the self-consistent soliton is employed as can be seen from 



®A test of the numerical treatment in the finite box is given in appendix A. 
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m/MeV cjrfig/MeV 2e^ai B^jMeY 

350 455 495 40 

400 368 421 53 

450 295 360 65 

500 225 303 78 

600 80 189 109 

TABLE II. The same as Tab. | in dependence on the constituent quark mass m for the 
self-consistent soliton and §2 = gi- 

Tab. J|. This feature is particularly astonishing because the valence quark energy de- 
creases simultaneously. 

These results for the diquark serve as input to the solution of the Faddeev equation 
( p. 37] ). In Fig. ^we display the resulting energy fig as a function of the scaling variable a 



for a constituent quark mass m = 450MeV. For comparison, the energy of three uncorre- 
lated valence quarks and the energy of a free diquark-quark pair are also shown. Again, in 



the entire region the diquark and quark are bound, although only weakly. In Tab. |T| we 
compare the quark-quark {Bgq = 2evai — ^diq) and diquark-quark {Bqd = e^ai + ^diq — ^b) 
binding energies to the total binding energy 

Btot = <i<^val —^B = Bqq + Bqd (4.1) 

for a diquark coupling constant g2 = Qi- We notice that the binding energy is carried to 
a large extent by the quark-quark correlations. While these have their maximal value for 
a vanishing soliton (a —>■ 0) the residual diquark-quark interaction develops a maximum 
for moderate soliton extension. More distinct is the situation (Tab. [TVD for a diquark 
coupling constant twice as large, i.e. g2 = 2gi. In particular, the total binding energy is 
strongly dominated by the diquark binding and is maximal for a vanishing soliton. 

In the past some authors have used instantaneous or static approximations to solve the 



integral equation [P^ ,P4 H26| , [42| . The instantaneous approximation corresponds to a quark 



exchange with vanishing relative energy = in eq. (|3.37|) . Then the different quark 



channels decouple and the Faddeev equation simplifies to a purely algebraic equation: 
r-o 1 •^^2 a 2I {val\Tdiq\K){K\Tdiq\val) 

[til - e^al - UJdiq] [{eval " UJdig - t^) - ^i^diq\ -^1^ — ^ = 0- (4-2) 

In the static approximation in addition the sum over the exchanged quark is restricted to 
the valence quark: 
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FIG. 3. Solution of the Faddeev equation (solid line), the energy of a free diquark-quark pair 
(short dashed line) and the energy of three uncorrelated valence quarks (long dash line) in dependence 
on the scaling factor a for a constituent quark mass m = 450MeV and g2 = Qi- 



a J^s/MeV Bgq/MeV Bqd/MeV Btot/MeV 

0.0 1221 175 6 181 

0.2 1222 173 6 179 

0.4 1215 171 8 179 

0.6 1010 136 66 202 

0.8 705 95 40 135 

1.0 452 65 24 89 

1.2 249 47 15 62 



TABLE III. The solution Jl^ of the Faddeev equation as well as the quark-quark {Bqg) and 
diquark-quark binding energy Bg^ and the total binding energy Btot of the three valence quarks for 
a constituent quark mass m = 450MeV and g2 = gi- 



ll 



a 




Bqq/MeY 


Bqd/MeV 


Btot/MeY 


0.0 


952 


433 


16 


449 


0.2 


957 


428 


17 


445 


0.4 


949 


401 


43 


444 


0.6 


787 


315 


110 


425 


0.8 


547 


169 


124 


293 


1.0 


374 


108 


59 


167 


1.2 


249 


75 


36 


111 



TABLE IV. The same as Tab. [Tj for the coupling constant §2 = 2gi. 

{val I Tdiq I val) {val \ Tdig \ val) _ ^ 

The presence of the energy denominator indicates that these approximations should be 
well suited for large quark and diquark masses. In fact, Tab. confirms this assessment. 
For a vanishing soliton the error is one per cent. When the soliton extension grows the 
valence quark as well as the diquark energy become smaller and therefore the error rises 
to about 20%. This emphasizes the importance of a complete solution of the Faddeev 
equation for a solitonic background. In addition, this table shows only a small difference 
between both approximations. In comparison to the complete solution of the Faddeev 
equation both treatments possess an enhanced binding energy. Hence we conclude that 
retardation effects in the full calculation have a repulsive character. Both approximations 
also elucidate the importance of the above described variable substitution, because the 
free solution iVL = ujdiq + e^ai is obvious for vanishing diquark-quark coupling {Tdiq 0). 
If we improperly had carried out the static approximation in equation (|3.36| ) we would 
have had two unphysical solutions ifl = 2ejjai and iQ = ±u!diq- 



[in — t^al — ^d: 



iq\ 



V. A HYBRID MODEL FOR BARYONS 

Given the results of the last section we now consider a soliton configuration with 
correlated valence quarks. For that purpose we replace the contribution ?>eyai of the three 
uncorrelated valence quarks in the energy functional ( |3.13D by the solution of the Faddeev 
equation. The total energy of this configuration is given by 

-^1^'' ~ ^sea + = Esol — Bfot < Esol{= Esea + '^^val) (5.1) 
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0.0 1221 1202 1201 

0.2 1222 1204 1202 

0.4 1215 1194 1189 

0.6 1010 935 924 

0.8 705 630 628 

1.0 452 394 393 

1.2 249 203 202 



TABLE V. The Solution of the Faddeev equation ( [3.37 ) as well as the solution of the static 



approximation Jl^*"*^ (|43| ) and the instantaneous approximation fi^"''*'* ( |4T2|) in dependence of the 
scaling factor a. 



1500 



1400 ' 
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>, 1300 
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1100 
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scaling factor a 

FIG. 4. The energy of the soliton with correlated valence quarks (solid line), the energy 
of the additive diquark-quark-soliton model (dash-dotted line) and the energy of the soliton with 
uncorrelated valence quarks (long dashed line) in dependence of the scaling factor a for a constituent 
quark mass 450MeV. 



19 



m (MeVj 




400 


450 


500 


600 


1 — t R /TV T T T\ 

^fo« (MeV) 


1057 


1140 


1142 


1108 


1011 


^sea (MeV) 


rv 







57U 


5o9 


606 


Ub (MeVj 


1057 


1140 


572 


539 


405 










0.90 


0.85 


0.80 


LOdiq (MeV) 


688 


725 


375 


356 


283 








0.56 


0.53 


0.47 


9A 


0.32 


0.32 


0.70 


0.68 


0.64 



TABLE VI. The total energy of the lowest energy configuration E^^i as well as the contribution 
from the Dirac sea Egea and from the correlated valence quarks il^ for §2 = gi- Furthermore 
displayed is the scaling parameter a, the effective energy LOdiq of the diquark, the isoscalar radius 
(rl^p)^/^ and the axial coupling constant qa- 



and is shown in Fig. ^ as a function of the scahng parameter a for the constituent quark 
mass of 450MeV. For comparison, the energy of the pure sohton ( p.l3| ) Egoi = Esea + '^^vai 
and the energy of the additive diquark-quark sohton model E^'^l'^^ = Esea + ^diq + ^vai 
defined in are also shown. As compared to the additive diquark-quark soliton picture 
the soliton with correlated valence quarks has a lower energy in the entire range shown. 
The minimum of this curve defines the hybrid model. Two other points in this figure 
also have a special meaning: First, the point a = on the same curve marks the model 
where a baryon is described as a bound state of a quark and a diquark with the meson 
fields fixed at their vacuum expectation values [^,^. Note, that this configuration only 
corresponds to a local minimum of the energy. The other important point is the minimum 
of Egoi at a = 1 which is nothing else than the self-consistent soliton solution without any 



valence quark correlations [|20|. Also the energy of this configuration is larger than the 
energy of the hybrid model. 

In Table we display results for the lowest energy configuration for different con- 
stituent quark masses. This table shows that for a large enough constituent quark mass 
the hybrid soliton in fact is the minimal energy configuration. The lower limit of stability 
for the hybrid soliton is around m ~ 430MeV. Below this value the energy of the baryon 
El^i = fi^ is carried completely by the valence quarks. Herein, the diquark possesses an 
effective energy of about TOOMeV. In comparison, for m ^ 450MeV the total energy of the 
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a 


„val 

9a 


„Fad 

9a 


diq 
9 A 


9a 


9A 


0.0 


0.325 


0.325 








0.325 


0.2 


0.325 


0.325 


0.001 


-0.001 


0.325 


0.4 


0.324 


0.323 


0.017 


-0.012 


0.328 


0.6 


0.273 


0.253 


0.478 


-0.033 


0.698 


0.8 


0.246 


0.239 


0.485 


-0.032 


0.692 


0.9 


0.240 


0.236 


0.483 


-0.018 


0.701 


1.0 


0.237 


0.234 


0.483 


0.005 


0.722 


1.1 


0.235 


0.233 


0.485 


0.035 


0.753 


1.2 


0.234 


0.232 


0.489 


0.074 


0.795 



TABLE VII. Contributions to the axial coupling constant in dependence of the scaling factor 
a. For comparison also the contribution of an uncorrelated valence quark is displayed. 



self-consistent soliton configuration is composed half by the valence quarks and half by 
the polarized sea. The effective diquark energy drops down to the range 280 — 380MeV. 

In Table we also present results for the isoscalar radius (rf^Q)^/^ and the axial 
coupling constant g^. These calculations]^ need some comments: The isoscalar radius is 
calculated as the spatial expectation value 

(r^^g) = J d^rr^p{r) = An J^^ drr^p{r). (5.2) 

This expression is only well-defined for localized baryon densities. For non-localized 
baryon densities the radius depends on the box size D which is the case when the soliton is 
absent[^. For the constituent quark mass m = 450MeV the isoscalar radius in the hybrid 
model (0.56fm) and the pure soliton model (0.54fm) are almost identical. We conclude 
from this slight increase that the sea is less attractive in the hybrid model despite of 
the stronger binding of the valence quarks. Note that both predictions are close to the 
experimental value (r|^Q)^/^ = 0.62fm. 

^In the appendix C we give a full description of the calculation of static observables in the 
hybrid model. 

^^Because of this volume divergence we state nothing in this case. 
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FIG. 5. The energy of the soliton with correlated valence quarks for different diquark coupling 
constants for a constituent quark mass 450MeV. 



More interesting are the results for the axial coupling constant 



Qa- Qa + 9a + 9a 



(5.3) 



consisting of contributions from the sea part Qa""^ the diquark g'^^'^ and the residual bound 
valence quark g^°''^ {cf. appendix C). These contributions are collected in Tab. [VI1| as a 



function of the scaling factor a. For comparison we also display the contribution of 
an uncorrelated valence quark. First, we observe only a small difference between a free 
and a correlated valence quark. As a function of the scaling factor we find two distinct 
regions: For small a (no soliton) we have g^"''^ ~ 0.33, for larger a (soliton region) the 
value decreases to g^"""^ ~ 0.24. On the other hand, the diquark contribution shows the 
opposite behavior. It vanishes in the case of a trivial background but it increases to more 
than twice the value of a single valence quark in the presence of the soliton. Finally, in 
both regions only small corrections occur, which originate from the sea contribution. The 
vanishing diquark contribution for a = also explains the small values for m = 350MeV 
and m = 400MeV in Tab. VI. In addition, these results indicate a missing component to 



the axial coupling constant: axial diquarks should also be included before comparing gA 
with the experimental value gA = 1-26 |^ |. 



When we increase the diquark coupling g2 we observe that stronger diquark corre- 
lations counteract the formation of a soliton. As can be seen from Fig. ^ the meson 
configuration, which minimizes E^^i , is shrinking with increasing g2. Starting at a critical 
coupling a path is open to a vanishing soliton. The displayed curves correspond to four 
special cases. The pure soliton is given for (72 = 0, whereas the second one {g2 = gi) 
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describes the Fierz-symmetric case. In the case g2 = 2gi the energy El^/ is identical to 
the physical nucleon mass [22] in the range where no soliton is present (a < 0.4). Finally, 
in the last case g2 = 3gi we have an additional symmetry - the Pauli-Giirsey symmetry 
H-l6| - giving the diquarks the same mass as the pions. However, this symmetry cannot 



be realized in nature since otherwise the baryon masses should be much smaller. 



VI. CONCLUSIONS AND OUTLOOK 



In this paper we have constructed a hybrid model for the description of baryons within 
the framework of the NJL model. Using functional integral techniques we have converted 
the pure quark model into an effective theory of mesons, diquarks and baryons. As a 
first step we have solved the Bethe-Salpeter equation for an S'-wave scalar diquark in the 
background of a soliton configuration. We have found that this diquark is kinematically 
stable against the decay into the two lightest quarks even in case the masses of the latter 
are strongly reduced by their interaction with the chiral soliton. In the next step we 
have integrated out the diquark fields, thereby inducing an interaction between quarks 
and the bound diquark through quark exchange. Using the ladder approximation for this 
exchange interaction a Faddeev type of equation has been derived, which includes the 
solitonic background field. This equation yields a three-quark bound state in the solitonic 
background. Using this result the hybrid has been constructed combining the three-quark 
bound state with the polarized sea. The energy of this hybrid is (for constituent quark 
masses m > 430MeV) smaller than either the energy of the soliton with three uncorrelated 
quarks or the energy of the three-quark bound state with no soliton present. Hence the 
formation of such a hybrid is dynamically preferred as compared to either picture alone. 
In addition, the soliton size of the hybrid is much smaller than the size of the soliton with 
uncorrelated quarks, and the Faddeev amplitude, i.e. the baryon wave function in terms 
of quarks, significantly deviates from the case without soliton. These results show that 
baryons are presumably very much like hybrids containing both, solitonic meson clouds 
and three-quark correlations, and that the implementation of both features is crucial for a 
proper description of baryons. Unfortunately, the question whether a baryon is dominated 
by the valence quark structure or the soliton is still an open question. E. g. if one includes 
for the soliton also (axial-) vector mesons the baryon current is completely carried by the 
polarized Dirac sea [^,0. To clarify the situation one should in addition include at the 
mean field level static diquark fields. The solution of such an improved soliton calculation 
would show whether there is still room for valence quark degrees of freedom beside vector 
meson degrees of freedom or whether they would exclude each other. 

Within the hybrid model we have calculated the isoscalar radius (rf^g)^/^ as well as 
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the axial coupling constant gA- While the isoscalar radius shows reasonable values in 
comparison with experiment the inclusion of axial diquarks seems to be important for 
the axial coupling constant. In order to calculate other observables one has to perform 
the last step in the construction of the hybrid model for baryons: one has to project the 
soliton solution to good spin and flavor states. 

In this work we have considered the NJL model which is well suited for low-energy 
properties of hadrons. Nevertheless, for energies higher than twice the constituent quark 
mass, threshold effects associated with the unphysical decay into free quarks disturb the 
calculation. Therefore, a generalization to bilocal chiral models seems necessary in order 
to effectively incorporate quark confinement. Within these models there can also exist 
confined diquarks |^9| on the one hand. On the other hand, soliton calculations have shown 
that inside a soliton propagating quarks [p0| -|5^ also can exist. Thus, the hybrid model 
introduced in this work presents a well suited starting point for further generalizations 
which appear necessary to deepen our understanding of the structure of baryons. 
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APPENDIX A: TEST OF THE NUMERICAL TREATMENT IN THE FINITE 

BOX 

In order to test the numerical treatment in the finite box we compare the solution of 
( ^.17] ) for a vanishing chiral angle 0(r) = with the solution in the continuum. The mass 
of the diquark is obtained from the Bethe-Salpeter equation [ p^ ] 



diq 



with 



B 



6momg2/gi 

' VP? + a(l — 







(Al) 



167r2 Jo 



167r2 



-1' A2 



A2 



(A2) 



When we compare the solutions of the Bethe-Salpeter equations in the continuum 
and in the spherical box ( |3.17| ) we observe a relatively large difference (c/. the 
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92/91 mdiq/MeY ujdiq/M-eV 



0.0 
0.5 
1.0 
1.5 
2.0 
2.5 
3.0 



900 
884 
736 
595 
463 
326 
140 



934 
904 
759 
625 
502 
380 
241 



934 
919 
778 
646 
527 
412 
288 



909 
893 
747 
608 
480 
350 
188 



TABLE VIII. The mass rudiq of a diquark in the continuum, the energy uj^iq of a diquark as 
the solution of the Bethe-Salpeter equation in a spherical box, w^^^ and w^^^ from the approximation 
of eq. ( [A^ ) and eq. ( |A4D respectively in dependence of 52/51- 



2"''^ and 3^'^ column in Tab. |VIII|) . However, we have to take into account that in 



a 



finite box the quarks always possess finite momenta and therefore an increased energy 



= ±-y/m^ + (nir/Dy. Now, we can consider two limits: First, in the limit of vanishing 
diquark coupling {g2 — ^ 0) a diquark consists of two free quarks {uJdiq = '^^vai)- From this 
we obtain the approximation 

= J^lq + K^/DY (A3) 



^diq — V '"'diq 

for the energy of a diquark in a spherical box. Second, we consider the limit of a large 
quark-quark coupling. In this case, the diquark is a highly correlated state which can be 
viewed as an independent particle, hence 



As can be seen from Tab. VIII and from Fig. ^ the diquark energy always lies within 



the two boundaries (^) and (^). In fact, in the limit of vanishing diquark coupling 

diq' 



{92 — ^ 0) the energy of the diquark Udiq is very well approximated by io^^ . The difference 



(2) 

at large coupling constants between uodiq and cu^j^ appears because of numerical effects at 
the boundary where the wave functions do not vanish completely. In the context of the 
hybrid model we consider bound states in the solitonic background. In this case, the wave 
functions are localized at the origin and the boundary effects are negligible. 
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1000 




FIG. 6. The energy of the diquark in a spherical box (solid line) in dependence of 52/51. the 

'£) and lower (wg) 



dashed lines corresponds to the upper (1^^]^ ) and lower (1^^^]) limit of the energy of the diquark 



APPENDIX B: THE BARYON PROPAGATOR 



The baryon propagator can be evaluated by comparison of eq. ( p.l9| ) and ( |2.12| ) 



(Bl) 



Using the definitions of ^ and ^ ( p.l3|) we can write 

'B )aP = 



/3 



(B2) 



y — A*G2iA^ A*G22A/3 J 



where we have used the inversion of the quark propagator 

Go 



G2I G22 



(B3) 



which is obtained from eq. (p. 14 ) by neglecting the baryon sources on the diagonal ele- 
ments. The elements of this matrix are given by 



Gn = (l + e?oA^~oA) 'e?o, 



G12 
G21 



6?oA l + 6?oA6?oA g^V, 



-1 ~ 



v-'g^ ^Kii + g^Ag^K) 'go, 



(B4) 



G22 = (1 + ^oA^^oA) ' goV. 
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APPENDIX C: CALCULATION OF STATIC OBSERVABLES 



The calculation of observables in the hybrid model separates in a valence quark part 
and a sea part: 



(CI) 



Because of the ansatz for the baryon wave function ( p.33|) 



{e^>0) 

\B) =Y^ ttaula, v) 



(C2) 



with 



(<:^>0) 



B{x,y) = {x,y\B) = / dudEY^ ^ a«^(E, , X4; y, 1/4) 



(C3) 



we can in addition separate the valence quark part in a diquark and a residual valence 
quark part: 



{0),ai = {B\0\B),,i 

a, (3 

(e^>0) ^ (e^,M>0) 

= {d)diq + {d)Fad. 

To take into account the internal quark structure of the diquarks we use the form 



(C4) 



= {d) = J VqVq q^Oqexp Q J qG \ 
= ^J VqVq qOqexp Q J qG^^q^ 
= ^ I VqVq exp Q I q (Q-^ + p6) q 



(C5) 



/3=0 



^Trlog(G-i+/36 



d(3 



f3=0 



where we have extended the operator O according to the Nambu-Gorkov formalism: 



O 



^ -foO ^ 



V 



-vdv^ 
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(C6) 



with O = VO'^'yoV~^. Expanding the expression (|C5| ) in powers of the diquark field the 
zeroth order term corresponds to (z/|0|/i) in the quark term whereas the second order 
term yields the diquark part (a|0|/5) in eq. (|C^) . 
The part of the residual valence quark is given by 



{0)Fad = I —Yl J2 a*^^{inB,'i'^)aaii{i^B,i^){j^\d\iJ,). 



The sea part is calculated as usual in the pure soliton model ^0 
{d)sea = ^(z/|6|z/)sgn(e^)erfc ( ^ ) • 



(C7) 



(a 



1. Baryon number B 



In the Nambu-Gorkov formalism the baryon number operator is defined as 

with B = Ic/Nc = lc/3. 



B 



^loB ^ 



V 



-V-foBV^ 



The term of second order in the diquark field is expressed as 



with the regularization function 



(C9) 



(CIO) 



(Cll) 



+ (e, - e^)^) .1 r (I, \ael + (1 - a)el + a{l - a)ujyA']) 
' daa — — 



4v^ 



ae^ + (1 — a)el + a{l — a)uj'^ 



5/2 



2. Baryon density 

For the calculation of the baryon densities we need the identity operator O = 1 as well 
as a spatial dependent variation variable (3 = (3{r) in eq. (|C5|) . The diquark contribution 
can be written as 

B,,^{t) = Y: Hir)^,ir)Rsitu; e., e,) (C12) 
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with the regularization function 



1 



A^/n Jo 



da ^ 



r (I, [ae^ + (1 +«(!-« VVA']) 



aef^ + (1 — a)ef. + a{l — a)c<j" 



3/2 



r (|, [ael + (1 - a)e2 + a(l - a)a7VA2] 
[ae2 + (1 - a)el + a(l - a)cu2]3/2 
[iuj + - e^] [-io; + - e^] ~ 

^ {lJ,\Tdig{tUj)\K}{K\Tdiq{-tUJ)W} 

r(|,[ae2 + (l-a)e2+a(l-«VVA2l 



(C13) 



ae^ + (1 — «)e^ + a{l — a)uj'' 



(/ilrdigFdjgl/i) 



3. Axial coupling constant 



Using the operator O 



.1^3 3 
3 



20 1 the axial coupling constant is given by 



(C14) 



with the regularization function i?3 defined in ( |C13|) . Note, that eq.( |C14]) is only valid 
for grand spin eigenstates 
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